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Abstract. The Fokas-Lenells (FL) equation arises as a model eqution which describes for nonlin- 
ear pulse propagation in optical fibers by retaining terms up to the next leading asymptotic order 
(in the leading asymptotic order the nonlinear Schrodinger (NLS) equation results). Here we present 
an explicit analytical representation for the rogue waves of the FL equation. This representation is 
constructed by deriving an appropriate Darboux transformation (DT) and utilizing a Taylor series ex- 
pansion of the associated breather solution, when certain higher-order nonlinear effects are considered, 
the propagation of rogue waves in optical fibers is given. 
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1. INTRODUCTION 

The Fokas-Lenells equation (FL)[H [2j d M 

kxt - iqxx + 2qx - q x qq* + iq = o, (l) 

is one of the important models from both mathematical and physical considerations. In eq.([T]), q 
represents complex field envelope and asterisk denotes complex conjugation, and subscript x (or t) 
denotes partial derivative with respect to x (or t) . The FL equation [TJ [3] is related to the nonlinear 
Schrodinger (NLS) equation in the same way that the Camassa-Holm equation is associated with the 
KdV equation. In optics, considering suitable higher order linear and nonlinear optical effects, the FL 
equation has been derived as a model to describe femtosecond pulse propagation through single mode 
optical silica fiber and several interesting solutions have been constrcuted [5]. 

Here, considering the physical significance of the FL equation and inspired by the importance of the 
recent interesting developments in the analysis of rogue waves of the NLS and the DNLS equations, 
we shall construct rogue wave solutions of the FL equation, by using the Darboux transformation 
(DT) [6], El [HI [9], [TUl [TT] . Our construction reveals that there exists a difference between the FL system 
and other integrable models, like the Ablowitz-Kaup-Newell-Segur (AKNS) |12[ [13] system and the 
Kaup-Newell (KN) [HJ E] system. 

Before discussing the DT of the FL system, let us briefly discuss the importance of rogue waves in 
mathematics and physics. Rogue waves have recently been the subject of intensive investigations in 
oceanography [151 [To] Tl\ [TBI R§\ [20] . where they occur due to modulation instability [21] [22[ [23[ [241 [25[ 
EH E7J [28j EH] , either through random initial condition [201 ED] or via other processes [T7 ] ITS1 ITS] \21 \ 
jMTM\M\M\M^^^M\\32\- The first order rogue wave usually takes the form of a single peak 
hump with two caves in a plane with a nonzero boundary. One of the possible generating mechanisms 
for rogue waves is through the creation of breatherspT| ITHl fT5] 1^ 1^ 1^ 1251 127^ 1551 RTTj which can be 
formed due to modulation instability. Then, larger rogue waves can emerge when two or more breathers 
collide with each other [33]. Rogue waves have also been observed in space plasmas [HI [34l ESJ [36], 
as well as in optics when propagating high power optical radiation through photonic crystal fibers 
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[3?1 l38j [39] . Considering higher order effects in the propagation of femtosecond pulses, rogue waves 
have been reported in the Hirota equation and the NLS-MB system [iQl UU 02]. 

This paper is organized as follows. Section 2, presents a simple approach to DT for the FL system, 
the determinant representation of the 1-fold DT, and formulae of anc | r PJ expressed in terms of 
eigenfunctions of the associated spectral problem. The reduction of DT to the FL equation is also 
discussed by choosing appropriate pairs of eigenvalues and eigenfunctions. Section 3 presents the 
construction of rogue wave solutions by using a Taylor series expansion about the breather solution 
generated by DT from a periodic seed solution with constant amplitude. Finally, in section 4 we 
conclude our results. 

2. Darboux transformation 

Let us start from the non-trivial flow of the FL (Fokas-Lenells) system [3] in the following modified 
form, 

iQxt - iQxx + 2q x - q x qr + iq = 0, (2) 
ir x t ~ ir xx - 2r x + r x rq + ir = 0, (3) 
which are exactly reduced to the FL eq.(pQ) for r = q* while the choice r = — q* would lead to eq.([I]) 
with change in sign of the nonlinear term. The Lax pairs corresponding to coupled FL eq.([2]) and ([3]) 
can be given by the FL spectral problem of the form [3] 

d x ip = ( JA 2 + QA)V = Ui/), (4) 
d t i; = (JA 2 +Q\ + V + V-iX- 1 + i JA- 2 )^ = V^, (5) 

with 

*-(*)■ '-(?!)■ 

n-i'-jr y-^{.i ft)- 

Here A, an arbitrary complex number, is called the eigenvalue(or isospectral parameter), and tp is 
called the eigenfunction associated with A of the FL system. Equations (|2|) and ([3]) are equivalent to 
the compatibility condition U t - V x + [U, V] = of (gj) and ©. 

Now, let us consider a matrix T of gauge transformation for the spectral problem ([!]) and ([S]) with 
the following form 

V> [1] =Tif>, (6) 

and then 

= U [1] U [1] = (T x + T U)T-\ (7) 

= V [1] V [1] = {T t + T V)T-\ (8) 

By cross differentiating ((7J) and ©, we obtain 

c/t^-yW, + [[/[!], y[ 1 ]] = r(c/ t -K + [c/,F])r- 1 . (9) 

This implies that, in order to prove eqs.fl2J) and ([5J are invariant under the transformation ©, it is 
crucial to construct a matrix T so that C/M and have the same forms as that of U and V. At the 
same time the old potentials (or seed solutions) (q, r) in spectral matrixes U and V are mapped into 
new potentials (or new solutions) (q^, r^) in terms of transformed spectral matrixes and V^. 
This newly obtained matrix T is a Darboux transformation of the FL system of eqs. (j^J) and (J3J). 

Based on the DT for the NLS[H \8\ and the DNLSP HOl ITT]. we suppose that a trial Darboux 
matrix T in eq.© is assumed to be in the form 

n-WM-(°S s ») + (V £)a-. do, 
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Here ai,di,bo,co are undetermined coefficients and they are function of (x, t), which will be expressed 
in terms of the eigenfunctions associated with Ai and A2 in the FL spectral problem and a_i and d-i 
are constants. For our further analysis, setting two eigenfunctions ipj as 



rj 



(11) 



After tedious but straightforward calculations, we get the one-fold Darboux transformation of the FL 
system as 

aiA + a_ 1 A _1 60 \ 



co 



Ti(A;Ai;A 2 ) 

and then the corresponding new solutions and rM are given by 

Jl] _ a -i , b o 



(12) 



a-i . wo [11 
a_i a_i 



a_i a_i 



(13) 



With the condition a_i = d_i = 1 and 
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Having obtained the values of all the unknown coefficients, we are now in a position to consider 
the reduction of the DT so that = (rM)*, then the DT of the FL equation is given. Under the 

h 



reduction condition q = r* with ip^ 



is an eigenfunction associated with the eigenvalue A&, 



then A; = X* k (k ^ I) and (pi = <f}-*, (fi = 4>k*- It is easy to check = (rM)* with the help of following 
choice of eigenvalues and eigenfunctions: 



Ai <-> tpi 



91 



, andAo 



AJ,o ^2 



<Pi 
4>\ 



(14) 



Thus, under the above choice of eigenfunctions and eigenvalues, the resulting T\ in eq. (|12p can be 
called as one-fold Darboux transformation of the FL equation. Further, taking one Pair of eigenvalues 
Ai = ot\ + i/3\ and A 2 = a\ — if3\ and their corresponding eigenfunctions, the final form of one-fold 
DT eq. (|12p reduces to a new solution of the form 



q + 



(A 2 2 -Ai 2 )0 m < 



A 2 Ai(-A 2 ^i</5i* + Ai0i0i*) 



(15) 



Notice that the denominator of has a nonzero imaginary part if f3± is not zero, which clearly 
indicates that the new solution gM { s a non-singular solution. 



3. ROGUE WAVES 

Using the results of DT constructed above, breather solutions of FL equation can be generated by 
assuming a periodic seed solution of eq. (|15|) . then we can construct the rogue waves of the FL equation 
from a Taylor series expansion of the breather solutions. 



For this purpose, assuming a and c as two real constants, then q = ce a> is a 

periodic solution of the FL equation, which will be used as a seed solution of the DT. Substituting 
the above form of q into the spectral problem eq.([4|) and eq.([5]), and using the method of separation 
of variables and the superposition principle, the eigenfunction ((pi,(pi) associated with Ai are given 
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explicitly. For simplicity, setting c = , taking ((f)i,<pi) back into eq. ([T5j) . then the 

a 

final form of the breather solution is obtained in the form 

mi f — 2aiPi(w\ cosh(ili) + W2sinh(Hi) + w^cos(H 2 ) + W4sin(H2)) 



qW = (- 



w§ cosh(//i) + Wf)Sinh(Hi) + u>7 cos(#2) + w%sin{H2) 



^a + 2a l *-2p l \ 2(j3 1 2 - ai 2 )t ^ 

) exp(i(ax + (a + 2 5 ))), (16) 

a ar 



where 



hi(—hs + 2y / /i3ai), ^2 = i\/ —hih 2 {—h^ + 2y / /i^ai), 



io 3 = 2qi/ii(-2qi + \/^3)) w 4 = — 2i\/— /ii/i2/3i(— 2ai + V^3 
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2ai/3i(ai 2 + /3i 2 )/ii(-2ai + \Zf^),w 6 = -2iaif3i{ai 2 + /?i 2 ) V-ZuM" 2 "! + v^), 



io 7 = /^(ax 2 + /3i 2 )/ii(-/i3 + 2^/feai), ws = -iai(a\ 2 + /3i 2 ) y / -/ii/i 2 (-/i3 + 2y / /i3~ai) 



i?l = — ^^((2a( ai 2 + A 2 ) 2 + ( ai 2 - ^ 2 ))t + 2a( ai 2 + £ 2 - 



h2 = a- 2(ai 2 + ft 2 ), /i 3 = a + 2(ai 2 - ft 2 ). 

From eq. (!16p . we can easily get the Ma breathers [26J (time periodic breather solution) and the Akhme- 
diev breathers |18[ [191 [25] (space periodic breather solution) solution. In general, the solution in 
eq. (!16p evolves periodically along the straight line with a certain angle between x and t axis. The 
dynamical evolution of |<?W| 2 in eq. (|16p is plotted in Figure 1. 
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Figure 1. Dynamical evolution of |gM | 2 (time periodic breather) with specific pa- 
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rameters a\ = 1, fi\ = — , a = The trajectory is a line at x = 0. 



Finally, a Taylor series expansion about the breather solution of the FL equation is used to construct 
the rogue wave of the FL equation in the following way. By letting a — > 2(a\ 2 + /3i 2 ) in ()16p . we obtain 
the rogue wave solution in the form 

(17) 



[1] s 2 

f ri l = ai — expsi 
S3 



with 



si = i(2xs + —t(As 2 - 2 13x2 + 3 +As + 1)), 
2s s 

s 2 = (t 2 (-/3i 2 + s)h + 16x 2 (-^ 2 + s)s 4 
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+ 8te(-/3i 2 + s)k 2 s + 2it(-/3i 2 + s)(12s 2 + l)s + 24ix(-^ 2 + s)s 3 - 3s 3 ), 
s 3 = (i 2 (-/3i 2 + s)h + 16x 2 (-/3i 2 + s)s 4 

+ 8te(-/3i 2 + s)/c 2 s - 2it(-/3i 2 + s)(4s 2 - l)s - 8ix(-/?i 2 + s)s 3 + s 3 )s, 
A;i = 64(-/3i 4 + s/3i 2 )(2/3! 4 - 2s/3i 2 + s 2 ) + 128/3i 8 - 256s/3i 6 

+ 192s 2 /3i 4 + (-64s 3 - 16s)/3i 2 + 8s 2 + 16s 4 + 1, 
k 2 = 12(-/?! 4 + s/3i 2 )s + 12s/3i 4 + (-12s 2 - 2)/3 x 2 + s + 4s 3 , 
s = ai 2 + /3i 2 . 

2 

By letting x — > oo, i — >■ oo, so k™| 2 — )• — ~ — . The maximum amplitude of |g™| 2 occurs 

(ai 2 + /3i 2 ) 2 

at t = and x = and is equal to - — ^— ,and the minimum amplitude of \q r w\ 2 occurs at 

(ai 2 +/3i 2 ) 2 

( _ Xa.' + ftV md :c _ T 3(12(a,» + ft')' + l) ^ „ ^ tQ „ Through ^ , ^ 

2^/3(ai 2 + 4/?i 2 )ai 8^/3(ai 2 + 4/?! 2 )ai 

Figure 3 of |glli| 2 , the main features (such as large amplitude and local property on (x-t) plane) of 
the rogue wave are clearly shown. 




Figure 2. Dynamical evolution of \q r 
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Figure 3. Contour plot of the wave amplitudes of k^l 2 with a\ 
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4. Conclusions 



Thus, in this paper, considering FL system of equation which describes nonlinear pulse propagation 
through single mode optical fiber, the determinant representation of the one-fold DT for the FL system 
is given eqs. (fl2j) and (fT3j) . By choosing pair of eigenvalues in the form A2 = A \ and assuming suitable 
seed solutions, the breather solution and rogue waves of the FL equation are derived in eq. (fl~6j) and eq. 
(|17p . Our results provide an alternative possibility to observe rogue waves in optical system. From the 
one- fold DT, it is interesting to observe that the DT of the FL system exhibits the following novelty 
in comparison with other integrable models like the AKNS and the KN systems: the DT matrix of 
the FL system has three different terms depending on A. 

Thus, the DT as well as the rogue wave of the FL system present novel features, in comparison 
with the DT and rogue wave solutions of the standard integrable systems like the AKNS and the KN 
systems. The construction of higher order rogue wave of the FL equation by using the determinant 
representation of the DT will be published elsewhere. 
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